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INTRODUCTION 


The  use  of  satellite  radar  altimeter  measurements  to  derive  an 
ocean  geoid  that  contains  geoid  undulation  information  at  much  shorter 
wavelengths  than  can  be  obtained  from  classical  analysis  of  satellite 
orbit  perturbations  has  been  postulated  for  a number  of  years,  e.g. 
references  1 and  2.  Some  experiments  in  this  direction  were  carried 
out  with  radar  altimeter  data  from  Skylab,  e.g.  references  3 and  4, 
but  the  amount  of  Skylab  radar  altimeter  data  and  the  lack  of  accuracy 
of  the  Skylab  orbital  ephemerides  precluded  geoid  determination  over 
any  appreciable  area. 

The  GEOS-3  satellite  was  launched  on  9 April  1975,  carrying  a radar 
altimeter  with  both  short  pulse  and  long  pulse  capability,  laser 
retroref 1 ectors  and  a Doppler  beacon  to  allow  precise  orbit  determination. 
GEOS-3  has  produced  radar  altimeter  data  of  both  good  quality  and  quantity 
for  over  a year.  Both  theoretical  studies,  reference  5 and  preliminary- 
evaluations,  reference  6,  of  the  GEOS-3  data  showed  that  the  short  pulse 
radar  altimeter  data  had  a lower  noise  level  than  the  long  pulse  data 
and  would  result  in  more  precise  derived  geofds. 

This  paper  develops  the  preliminary  derivation  of  a local  geoid 
in  the  Eastern  Indian  Ocean  from  94  passes  of  short  pulse  GEOS-3  radar 
altimeter  data  received  by  the  DoD  telemetry  station  at  Perth,  Australia. 
This  data  was  taken  over  a 34  day  period  in  July-August  1975.  The 
orbital  ephemerides  used  in  the  analysis  were  obtained  from  Doppler 
tracking  data  and  have  an  estimated  accuracy  of  2 meters  in  the  radial 
component,  reference  7. 

This  preliminary  derivation  of  the  geoid  does  not  consider  tide 
effects  on  the  altimeter  measurements  which  could  introduce  uncertainties 
of  about  0.5  meters.  Also,  the  standard  tropospheric  range  correction 
used  on  the  data  could  result  in  errors  in  individual  measurements  of 
up  to  25  centimeters.  A technique  for  reducing  all  the  biases  by 
analyzing  the  radar  altimeter  measurements  at  track  intersections  is 
employed.  A later  analysis  is  planned  that  will  include  data  from 
other  telemetry  stations  in  the  area  whose  coverage  overlaps  the  Perth 
coverage  and  extends  the  total  time  span  of  the  data  set  to  several 
months.  It  will  include  tidal  corrections  to  the  data  before  the  bias 
removal  and  geoid  derivation  is  carried  out. 

The  preliminary  GEOS-3  geoid  derived  in  the  paper  is  compared  with 
the  Marsh  and  Chang  1976  Geoid,  ref.  8.  Correlations  of  larger 
features  in  the  GEOS-3  geoid  with  bottom  topography  are  also  given. 


PROCEDURE 


Averaging  at  Intersections 

Ninety-four  passes  of  short  pulse  radar  altimetry  data  were  taken 
by  the  Perth,  Australia  telemetry  station  and  processed  at  NSWC/DL 
using  precise  short-arc  Doppler  derived  satellite  ephemerides.  For 
each  pass  time,  latitude,  longitude,  geoid  height  and  along-track 
deflections  of  the  vertical  were  computed  at  0.1  second  time 
intervals.  Details  of  this  processing,  including  the  Wiener  filter 
used,  may  be  found  in  reference  9.  Nine  of  the  passes  contained  data 
too  fragmented  for  efficient  operation  of  the  Wiener  filter  and  were 
eliminated  from  the  following  analysis.  The  remaining  85  passes  were 
used.  The  ground  tracks  of  these  passes  are  given  in  Figure  1.  The 
differences  in  geoid  height  at  each  of  the  716  intersections  of 
these  passes  are  given  in  Table  1.  The  average  difference  of  geoid 
height  over  all  the  intersections  is  75  cm  with  a standard  deviation 
with  respect  to  the  mean  of  2.3  meters.  The  maximum  difference  is 
7.83  meters.  As  a first  approximation  to  the  geoid,  the  average  of 
the  two  geoid  heights  at  each  intersection  is  taken  as  the  geoid 
height  and  contours  of  constant  geoid  height  plotted.  Since  along- 
track  deflection  of  the  vertical  is  available  for  each  track  at  each 
intersection,  the  surface  gradient  at  each  intersection  may  be 
computed.  This  is  used  as  an  interpolation  aid  in  plotting  the  geoid 
height  contours.  A geoid  derived  in  this  manner  is  given  in 
Figure  2.  The  +'s  in  Figure  2 represent  the  location  of  the 
intersections.  The  contour  interval  in  Figure  2 is  2 m.  The  General 
Purpose  Contouring  Program  by  California  Computer  Products  was  used 
to  generate  the  contours  in  Figure  2 and  all  subseguent  contour  plots. 

Bi as  Remova 1 

If  one  studies  the  intersections  in  Table  1,  it  is  seen  that  some 
passes  have  all  positive  differences,  e.g.  descending  pass  1580,  or 
all  negative  differences,  e.g.  descending  pass  1837.  Since  pass  1580 
is  intersected  by  19  other  passes  taken  over  a time  span  of  33  days 
and  pass  1837  is  intersected  by  16  other  passes  taken  over  a time 
span  of  32  days,  one  could  assume  that  different  tidal  phases, 
different  deviations  from  the  standard  troposphere,  different  orbit 
uncertainties  on  the  crossing  passes  would  be  sampled.  Therefore 
the  intersections  would  be  expected  to  be  more  random  unless  the  single 
pass,  i.e.  1580  or  1837,  was  biased  with  respect  to  the  other  passes. 
Using  this  assumption  each  pass  in  Table  1 has  been  adjusted  so  that 
the  sum  of  the  intersection  differences  with  the  crossing  passes  have 
been  minimized.  The  mean  bias  adjustment  is  0 cm;  thus  the  mean 
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geoid  is  not  changed.  The  standard  deviation  of  the  biases  is  2.2  m. 
The  amount  of  the  adjustment  of  each  pass  is  given  in  Table  2.  The 
intersection  differences  after  applying  the  adjustments  are  presented 
in  Table  3,  where  the  average  difference  summed  over  all 
intersections  is  now  9 cm  with  a standard  deivation  of  83  cm  with 
respect  to  the  mean  and  a maximum  difference  of  3.15  m.  The  geoid 
heights  are  again  taken  as  the  midpoint  at  each  intersection  and  a 
contour  map  is  made  as  before  using  the  surface  gradient  at  each 
intersection  to  aid  in  interpolation. 


RESULTS 


Compari son  With  the  Marsh  and  Chang  1976  Geoid 

The  preliminary  eastern  Indian  Ocean  geoid  derived  from  the 
GEOS-3  measurements  using  the  bias  removal  technique  described 
above  is  given  in  Figure  3.  The  +'s  plotted  on  Figure  3 indicate  the 
locations  of  the  intersections.  Figure  4 is  a plot  of  the  Marsh  and 
Chang  1976  Geoid.  A contour  plot  of  the  differences  between  the 
GEOS-3  Geoid  after  bias  removal  and  the  Marsh  and  Chang  1976  Geoid 
is  given  in  Figure  5.  The  contour  interval  in  Figures  3,4,  and  5 is 
2 n.  In  general  the  agreement  is  good,  with  major  features  such  as 
the  through  running  from  the  southwestern  tip  of  Australia  toward  the 
northwest,  and  the  saddle  point  south  of  Australia  comparing  well. 

The  maximum  difference  between  the  two  geoids  is  7 m.  To  further 
compare  the  two  geoids, the  along-track  geoid  heights  (after  bias 
removal)  from  pass  1419  and  descending  pass  1724  are  compared  in 
Figures  6 and  7,  with  the  Marsh  and  Chang  1976  geoid  heights  along 
these  tracks.  Again  the  agreement  is  quite  good,  with  the  GEOS-3 
along-track  data  showing  more  fine  structure  than  was  apparent  in  the 
o/erall  geoid  plots,  since  the  GEOS-3  contour  plots  used  information 
only  at  the  intersections  instead  of  at  every  point  along-track. 

Both  the  GEOS-3  Geoids  and  the  Marsh  and  Chang  1976  Geoid  were 
referenced  to  an  ellipsoid  with  a semi-major  axis  of  6378.135  km 
and  a flattening  reciprocal  of  298.26. 

Compari son  with  Bottom  Topography 

Comparison  of  the  GEOS-3  Geoid  was  made  with  larger  bottom  topography 
features  of  the  eastern  Indian  Ocean.  Figure  8 is  a contour  map  of  a 
small  area  of  the  geoid  with  a one  meter  interval  between  contour  lines. 
Figure  9 is  a three  dimensional  projection  of  the  same  area  shown  in 
Figure  8.  Comparison  of  Figures  3,  8 and  9 with  ocean  bottom  charts 
shows  good  agreement  with  large  features.  An  attempt  was  made  to 
correlate  the  GEOS-3  Geoid  with  finer  scale  bottom  topographic  features 
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shown  on  a deta i 1 ed  charts  supplied  by  the  Naval  Oceanographic  Office 
without  much  success,  probably  due  to  the  smoothing  of  the  GEOS-3 
Geoid  by  using  data  only  at  intersections  in  its  construction.  Much 
success  was  had  by  making  along- track  comparisons  where  the  short 
wavelength  information  is  not  lost.  Figure  10  shows  an  along-track 
comparison  of  bottom  topography  and  geoid  heights  for  pass  174G. 

Geoid  heights  are  plotted  for  the  GEOS-3  Geoid  after  bias  removal  as 
well  as  the  along-track  satellite  data.  The  lack  of  short  wavelength 
information  in  the  GEOS-3  Geoid  is  evident.  Correlation  of  the  three 
curves  is  very  good.  The  sea  mount  near  the  middle  of  the  pass  shows 
up  clearly  in  the  two  geoid  height  curves.  The  sea  mount  is  located 
approximately  at  lattitude  -2232  and  longitude  10436.  The  geoid 
undulations  caused  by  this  sea  mo'mt  can  be  seen  in  figures  2, 3, 4, 8, 
and 


CONCLUSIONS 


Careful  processing  of  GEOS-3  radar  altimetry  data  with  good  orbit 
ephemerides  can  result  in  a derived  ocean  geoid  that  is  accurate  to 
about  one  meter.  It  is  expected  that  greater  accuracy  will  be 
obtained  in  the  GEOS-3  Geoid  with  the  analysis  of  additional  data  in 
the  area  from  other  telemetry  stations  whose  coverage  partially 
overlaps  Perth  and  which  will  provide  data  over  a longer  time  frame. 
Also  the  inclusion  of  tide  and  sea  state  corrections  and  a bias  rate 
term  in  the  analysis  should  result  in  greater  accuracy  in  the  derived 
geoid.  The  use  of  along-track  data  (instead  of  using  only  intersection 
data)  should  improve  the  accuracy  and  the  short  wavelength  structure  of 
the  GEOS-3  Geoid. 
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3.73 

2.14 

1.93 

2.87 

3. 84 

. --  — ■ 


TABLE  2 


BIASES  (METERS) 

Ascending  Passes  Descending  Passes 


Pass  No. 

Bias 

1347 

2.0255 

1349 

1.4120 

1362 

.8006 

1376 

2.6962 

1390 

1.1268 

1391 

-2.6197 

1404 

.1451 

1405 

-1.5824 

1406 

2.7338 

1418 

-2.8210 

1419 

-1.1408 

1420 

1.5730 

1433 

2.9797 

1434 

1.5557 

1447 

2.3771 

1561 

2.0308 

1575 

3.6217 

1633 

.1213 

1646 

1.9914 

1647 

-.0689 

1660 

-.8489 

1661 

-.6609 

1662 

.3416 

1674 

-3.9360 

1675 

-1.5822 

1676 

.7676 

1689 

-.0942 

1703 

2.9032 

1704 

1.0856 

1717 

-.6313 

1718 

-.1651 

1731 

-4.4150 

1732 

-1.6223 

1745 

2.5125 

1746 

-.8584 

1747 

2.3044 

1759 

-3.2807 

1775 

.7127 

1803 

-.5055 

1831 

2.3792 

Pass  No. 

Bias 

1354 

2.0087 

1368 

2.8957 

1382 

3.6436 

1383 

-.7709 

1396 

.7984 

1397 

-.8719 

1410 

-.9172 

1411 

-.5216 

1424 

-2.4272 

1425 

1.5748 

1426 

-.7826 

1566 

-2.7478 

1568 

-1 .2155 

1580 

-2.8077 

1581 

1 .4333 

1624 

1 .1353 

1637 

-3.2068 

1651 

-2.2751 

1652 

1.2811 

1653 

-1.3553 

1665 

-.7464 

1666 

-.2933 

1667 

-.7286 

1680 

-2.3492 

1681 

-.2807 

1682 

-2.2569 

1694 

-2.5601 

1696 

.0629 

1708 

-2.4820 

1709 

3.1803 

1710 

-.9380 

1722 

-1 .3373 

1723 

1.0554 

1724 

-2.2592 

1736 

-1 .6476 

1737 

-1.1118 

1738 

-.5532 

1739 

-2.2550 

1752 

4.2271 

1809 

.5532 

1837 

4.4850 

17 


DISTRIBUTION 


r— 


Defense  Documentation  Center 
Cameron  Station 

Alexandria,  Virginia  22314  (12) 

Di rector 

Defense  Mapping  Agency 
U.S.  Naval  Observatory 
Bldg.  56 

Washington,  D.C.  20360 
Attn:  0.  W.  Williams 

Charles  Martin 
P.  M.  Schwimmer 
D.  A.  Rekenthaler 

Naval  Oceanographic  Office 
Washington,  D.C.  20390 
Attn:  T.  Davis 

J.  Hankins 

Pentagon 
OAD/SNAS 
Room  3E139 

Washington,  D.C.  20301 
Attn:  Leroy  Kuykendahl 

Pentagon 

JESCG 

Room  3E139 

Washington,  D.C.  20301 
Attn:  Capt.  H.  Bixby 

Office  of  the  Oceanographer 
Hoffman  2 
200  Stoval  Street 
Alexandria,  Virginia  22332 
Attn:  Capt.  D.  Brown 

Di rector 

Naval  Research  Laboratory 
Washington,  D.C.  20360 
Attn:  B.  S.  Yaplee 

V.  E.  Noble 
A.  Shaprio 


19 


Oceanographer  of  the  Navy 

Hoffman  2 

200  Stovall  St. 

Alexandria,  Virginia  22332 

Strategic  Systems  Project  Office 
Dept,  of  Navy 
Washington,  D.C.  20390 
Attn:  SP20 

SP23 
SP24 

Defense  Mapping  Agency 
Aerospace  Center 
St.  Louis,  Missouri  63118 
Attn:  R.  Ballew 

0.  Finkland,  Code  GDG 

Defense  Mapping  Agency 
Topographic  Center 
Washington,  D.C.  20315 
Attn:  R.  Smith 

Defense  Mapping  Agency 
Hydrographic  Center 
Washington,  D.C.  20360 
Attn:  Maj.  J.  Jerome 

Air  Force  Geophysical  Research  Laboratory 
L.  G.  Hanson  Field 
Bedford,  Mass.  01730 

Attn:  George  Hadgigeorge 

U.S.  Army  Topographic  Laboratory 
Ft.  Belvoir,  Virginia  22060 
Attn:  Armando  Mancini 


NASA 

Wal lops  FI ight  Center 
Wallops  Island,  Virginia  23337 
Attn:  J.  T.  McGoogan 

H.  R.  Stanley 
W.  F.  Townsend 


NASA 

Goddard  Space  Flight  Center 
Greenbelt,  Md.  20771 
Attn:  J.  W.  Siry 

F.  0.  Vonbun 
D.  E.  Smith 

J.  G.  Marsh  (2) 

NASA 

600  Independence  Ave.  , S.W. 

Washington,  D.C.  20546 
Attn:  F.  Williams 

S.  W.  McCandless 

NOAA-National  Ocean  Survey 
6001  Executive  Blvd. 

Rockville,  Maryland  20852 
Attn:  B.  H.  Chovitz 

B.  Douglas 

NOAA 

Pacific  Marine  Environmental  Lab. 

Seattle,  Washington  98105 
Attn:  John  R.  Apel 

NOAA-AOML 

15  Rickenbacker  Causeway 
Miami,  Florida  33149 

Attn:  H.  Michael  Byrne 

NOAA 

NESS/SPOC  Group 

5100  Auth  Road 

Camp  Springs,  Md.  20031 

Attn:  John  W.  Sherman  III 

NOAA-ERL 

Boulder,  Colorado  80303 
Attn:  Leonard  Fedor 

Battel le  Columbus  Laboratories 
505  King  Avenue 
Columbus,  Ohio  43201 
Attn:  A.  G.  Mourad 


21 


Jet  Propulsion  Laboratory 
4800  Oak  Grove  Drive 
Pasadena,  California  91103 
Attn:  James  A.  Dunne 

Alden  Loomis 

Weapons  Research  Laboratory 
Box  2151  GPO  Adelaide 
South  Australia  5001 
Australia 

Attn:  Garth  A.  Morgan 

Scripps  Institution  of  Oceanography 
LaJolla,  California  92093 
Attn:  Robert  Stewart 

University  of  New  South  Wales 
P.0.  Box  #1 

Kensington  NSW  2033,  Australia 
Attn:  R.  S.  Mather 

Smithsonian  Astrophysical  Observatory 
60  Garden  Street 
Cambridge,  Mass.  C2138 
Attn:  E.  M.  Gaposchkin 

City  University  of  New  York 

institude  of  Marine  and  Atmospheric  Sciences 

675W.  252  Street 

Bronx,  N.  Y.  10471 

Attn:  Prof.  Willard  J.  Pierson 

Woods  Hole  Oceanographic  Institution 
Woods  Hole,  Mass. 

Attn:  Dr.  John  Whitehead 


Local : 

DK-10 

(50) 

DK-50 

(10) 

DK-01 

DK-05 

DX-21 

(2) 

DX-222 

(6) 

DX-40 

22 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  of  This  PAGE  (When  Dmta  Entered) 


REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 

1 REPORT  NUMBER  

TR-3668  k 

2.  GOVT  ACCESSION  NO. 

/ ? 

3.  fl^ClP  1 E N T ’ S CATALOG  NUMBER 

4 TITLE  (and  Subtitle) 

Preliminary  Eastern  Indian  Ocean  Geoid  From  GEOS-3 

5.  TYPE  OF  REPORT  A PERIOD  COVERED 

Data  , 

i 7 . - — 

6 PERFORMING  ORG.  REPORT  NUMBER 

7 AU  T HO R( a) 

Samuel  L. /Smi th , III  / 
Alan  C...  Chappell  / 

8 CONTRACT  OR  GRANT  NUMBERfaJ 

9 PEffFORMING  ORGANIZATION  NAME  AND  ADDRESS 

Naval  Surface  Weapons  Center 
Dahlgren  Laboratory 
Dahlgren,  Virginia  22448 

10.  PROGRAM  ELEMENT.  PROJECT,  TASK 
AREA  A WORK  UNIT  NUMBERS 

DMA  PF  63701 B/3204/240 

II.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

Defense  Mapping  Agency,  ATTN:  STT  ( / -■ 

Bldg.  56,  U.S.  Naval  Observatory 

"achinntnn  HP  ? H ”3  f)  R 


4 MONITORING  AGENCY  NAME  A ADDRESS (ti  different  from  Controlling  Office)  | 15  SECU 


16  DISTRIBUTION  STATEMENT  (of  thle  Report) 


15a.  DECLASSIFICATION  DOWNGRADING 
SCHEDULE 


Approved  for  public  release;  distribution  unlimited. 


17  DISTRIBUTION  STATEMENT  (of  the  ebetract  entered  In  Block  20,  If  different  from  Report) 


19  KFY  WORDS  (Continue  on  reverse  aide  If  necessary  and  Identify  by  block  number) 


Geoid,  Indian  Ocean,  GEOS-3,  Satellite  Radar  Altimetry,  Bottom  Topography 


20  A^S1  RACT  (Continue  on  reverse  side  If  necessary  and  Identify  by  block  number) 

''Ninety-four  passes  of  GEOS-3  short  pulse  radar  altimetry  data  taken  in 
July-August  1975  over  the  eastern  Indian  Ocean  (by  the  DoD  telemetry  station 
at  Perth,  Australia)  have  been  analyzed  to  arrive  at  a preliminary  ocean  geoid. 
The  self  consistency  of  the  data  at  track  intersections  has  a mean  geoid  height 
difference  of  75  cm  with  a standard  deviation  of  2.3  meters  with  normal  pro- 
cessing. Application  of  bias  removal  techniques  reduced  the  mean  geoid  height 
difference  to  9 cm  and  improves  the  self  consistency  of  the  results  to  a 


DD  1473  EDITION  OF  I NOV  S5  IS  OBSOLETE 

> j*n  ti  mr  AQCTCicn 


S N 0102-  LF*  014-  6601 

_ ^ 


secunity  classification  of  this  page  twhtn  n«r»  imnM) 


standard  deviation  of  less  than  1 m.  A comparison  of  the  GEOS-3  Geoid  is  made 
vith  the  NASA  Marsh  and  Chang  1976  Geoid.  The  GEOS-3  Geoid  correlates  well 
Aiith  the  larger  features  of  the  local  bottom  topography. 


* 


